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Abstract — This article^ studies disruption tolerant networks 
(DTNs) where each node knows the probabilistic distribution of 
contacts with other nodes. It proposes a framework that allows 
one to formalize the behaviour of such a network. It generalizes 
extreme cases that have been studied before where (a) either 
nodes only know their contact frequency with each other or (b) 
they have a perfect knowledge of who meets who and when. 
This paper then gives an example of how this framework can be 
used; it shows how one can find a packet forwarding algorithm 
optimized to meet the 'delay/bandwidth consumption' trade-off: 
packets are duplicated so as to (statistically) guarantee a given 
delay or delivery probability, but not too much so as to reduce 
the bandwidth, energy, and memory consumption. 

I. Introduction 

Disruption (or Delay) Tolerant Networks (DTNs, [1]) have 
been the subject of much research activity in the last few 
years, pushing further the concept of Ad Hoc networks. 
Like Ad Hoc networks, DTNs are infrastructureless, thus the 
packets are relayed from one node to the next until they 
reach their destination. Moreover, in DTNs node clusters can 
be completely disconnected from the rest of the network. In 
this case, nodes must buffer the packets and wait until node 
mobility changes the network's topology, allowing the packets 
to be finally delivered. 

A network of Bluetooth-enabled PDAs, a village intermit- 
tently connected via low Earth orbiting satellites, or even an 
interplanetary Internet ([2]) are examples of disruption tolerant 
networks. 

The atomic data unit is a group of packets to be delivered 
together. In DTN parlance, it is called a message or a bundle; 
we use the latter in the following. 

Routing in such networks is particularly challenging since 
it requires to take into account the uncertainty of mobiles 
movements. The first methods that have been proposed in the 
literature are pretty radical and propose to forward bundles in 
an "epidemic" way ([3], [4], [5]), i.e., to copy them each time 
a new node is encountered. This method of course results in 
optimum delays and delivery probabilities, at the expense of an 
extremely high consumption of bandwidth (and, thus, energy) 
and memory. To mitigate those shortcomings, the epidemic 
routing has been enhanced using heuristics that allow the 
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propagation of bundles to a subset of all the nodes ([6], [7], 
[8]). 

Since node's buffer memory is not unlimited, a cache 
mechanism has been proposed, where the most interesting 
bundles are kept (i.e. those that are likely to reach their 
destination soon) and the others are discarded when the cache 
is full ([9], [10], [11], [12], [13], [14]). Those schemes must 
thus guess when a bundle will reach its destination, which 
is most of the time computed thanks to frequency contact 
estimation (which reflects the probability that two given nodes 
meet in the future). 

Few papers explore how the expected delay could be more 
precisely estimated (notable exceptions are [15], [16]). It has 
been proved ([17]) that a perfect knowledge of the future node 
meetings allows the computation of an optimal bundle routing. 

This short overview emphasizes two shortcomings: 

> Certain networks might be highly predictable (e.g. nodes 
are satellites and links appear and vanish as they revolve 
around their planet), others are much more chaotic. Previ- 
ous work suppose either that nodes contacts are perfectly 
deterministic and known in advance, or that only the 
contact frequency is known for each pair of nodes. We 
propose to generalise these approaches and suppose that 
each node knows a probability distribution of contacts in 
the (near) future. 

• [5] underlines the tradeoff between bundle delivery guar- 
antees and bandwidth/energy consumption: copying the 
bundles is costly since, in mobile networks, those re- 
sources are both scarce. Current schemes use a cache 
mechanism that ensures each node only receives the most 
relevant bundles, which somehow mitigates this problem, 
but does not provide any rationale, except the need to 
cope with mobiles limited memory. We propose to route 
the bundles according to the delivery or delay guarantees 
required by the user, thus only duplicating packets when 
it is beneficial. 

This paper is organised as follows. Section [TT] presents a 
way to model the contacts between the nodes of a predictable 
network. Sections [XT] and [IV] show how the end-to-end delay 
of bundles can be predicted. Sections [V] and [Vl] give a routing 
algorithm that allows to deliver bundles in a manner that meets 
a given guarantee. Section [VTI1 concludes. 
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Fig. 1. Contact profile of a node pair over a month: example. The height of a bar gives the probability that two nodes meet (at least once) during the 
corresponding 12-hour time period. Here, nodes are supposed to meet at the beginning of each week, but the exact day is unknown. 
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Fig. 2. First contact probability distribution corresponding to the contact profile figure [T] (Each bar corresponds to a 12-hour period.) 



II. Predictable future contacts 

The network is composed of a finite set of wireless nodes Af 
that can move and thus, from time to time, come into contact. 

In the sequel, a contact between two nodes happen when 
those nodes have setup a bi-directional wireless link between 
them. A contact is always considered longenough to allow all 
the required data exchanges to take placaj. 

A. Contact profiles 

We expect the mobiles motion to be, to a certain extent, 
predictable, yet obviously the degree of predictability varies 
from one network to another. Sometimes nodes motion is 
known in advance because they must stick to a given schedule 
(e.g. a network of buses) or because their trajectory can easily 
be modelled (e.g. nodes embedded in a satellite). Other net- 
works are less predictable, yet not totally random: colleagues 
could be pretty sure to meet every day during working hours, 
without any other time guarantee. Mobile nodes behaviour 
could also be learnt automatically so as to extract cyclical 
contact patterns. 

We therefore suppose that each node pair {a,b} C Af can 
estimate its contact probability for each time step in the near 
future. We call it a contact profile and denote it C a b ■ N — > 
[0,1]. The time step duration should be chosen small compared 
to the expected network's end-to-end delay. Figure Q] gives an 

2 This is a major difference with [17] which does not neglect bundle 
transmission times. 



hypothetical contact profile. In the following, we suppose the 
profile known for each node pair. 

Contact profiles can easily represent situations usually de- 
picted in the literature: 

• A constant profile C a b{t) — k describes a node pair 
that only knows its contact frequency. For example, the 
profile C a b(t) — 1/30 (probability of contact per day) 
corresponds to two nodes a and b meeting once a month 
on average. 

« Perfect knowledge of nodes meeting times results in a 
profile made of peaks: Vt S N : C a b(t) G {0, 1}. 
In practice, unknown contact profiles could be replaced by 
a null function to get a defensive approximation of their 
behaviour. 

The following sections aim at studying how bundles prop- 
agate from one node to another in a network whose nodes' 
contact profiles are known. 

B. First contact distribution 

It is easy to deduce the probability distribution of a (first) 
contact at time t between nodes a and b G Af given their profile 
C a b', we denote this distribution d ab - Since the probability of a 
first contact at time t is the probability of meeting at time step 
t times the probability not to meet at time steps 0, 1, . . . ,t— 1, 
we have: 

t-i 

d ab (t) = C ab (t) U(l-C ab {i)) Va,6 G Af, Vt G N (1) 

i=0 
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Fig. 3. The contact probability density -D t,(9, •) matching the contact profile given in figure[T] 



The distributions domain is N since contact profiles have 
been defined using discrete time steps. We extend the distri- 
butions to R to get rid of this artifact. Notice that d a b is not a 
well-defined probability distribution since its integral over its 
domain is not equal to 1: two nodes might never meet. Those 
considerations directly lead to the definition of the first contact 
distribution set. 

Definition 1: The first contact distribution set, C, is the set 
of function^ / : R+ -> R+ such that J °° f(x) dx < 1. 

Contact profiles have a shortcoming: they do not allow us 
to express contact interdependencies; for example, they cannot 
model that two nodes are certain to meet during the weekend 
without knowing exactly which day. First contact distributions 
have no such limitations. Therefore, when it is possible, one 
could find preferable to generate them directly without relying 
on contact profiles. 

Figure |2] gives the d a b distribution corresponding to the 
contact profile C a b depicted in figure Q] 

Notice that if a bundle is delivered directly from a to b, 
knowing the first contact distribution allows an easy verifica- 
tion of a large spectrum of guarantees, such as the average 
delay or the probability of delivery before a certain date. 

III. Delivery distributions 

A. Definition 

First contact distributions can be generalized to take into 
account the knowledge that no contact were made before a 
certain date. 

Let D a b(T,t) be the probability distribution that a and b 
require a delay of t time steps to meet for the first time after 
time step T. Since these distributions will be the building 
blocks that allow us to compute when a bundle can be 
delivered to its destination, we call them delivery distributions. 
D a b can directly be derived from the contact profile C a b- 

T+t-1 

Dab(T,t)=Cab{T + t) Y[ (1-C ab (l)) 
i=T 

Vo, b e AT, VT, t G N (2) 

3 R~'~ denotes the set of positive reals. 




Fig. 5. The D a ),(T,t) function matching the contact profile given in 
figure Q] 
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Definition 2: The delivery 
holds all the functions / 
VT G M+ : r o °° f(T, x) dx < 1. 
Notice the inequality. 

The D a b(T,t) distribution corresponding to the contact 
profile given in figure [T] is plotted in figure [5] Figure [3] plots 
the function D a b(9, •) (i.e. a section of D a b(T, t) in the T = 9 
plane); the D(T, ■) functions of course belong to C (VT > 0). 

Notice that D a b(T, •) is the expected delivery delay distribu- 
tion for a bundle sent directly from a source a to a destination 
6 if a decides to send it at time T. 

B. Order relation on distributions 

We define an order relation between first contact distri- 
butions. Intuitively, this relation allows us to compare two 
distributions to find which one represents more frequent or 
predictable contacts. A rigorous definition is given below. 

Definition 3: The first contact distributions di G C is 
greater (or equal) than di £ C (denoted di >z cfo) if and 
only if: 



\/x>0: / d x {t)dt > / d 2 (t)dt 
Jo Jo 



(3) 



This relation is a partial order (but not a total order as 
there exist d\, di G C such that neither d\ b $2 nor d\ < c?2)- 
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Definition [3] specifies when two 
contact distributions d\ , d 2 are such 
that d\ >z d,2- The plots show 
two distribution examples (left- 
hand plot) and their cumulative 
function (denoted d* and d|, right- 
hand plot). We have d\ > d 2 iff 
Mt > : df(t) > d* 2 (t). Here, 
neither d\ >z d 2 nor d2 cz di hold. 

The distribution sup{di, d 2 } (left-hand plot) is called supremum (or least 
upper bound). Its cumulative function is the maximum of the d* and d 2 
functions; its distribution is the derivative of the cumulative function. 
By definition of sup{d\, d 2 }, if d >; d\ and d >z d 2 , then 
d >z sup{di, d 2 } (Vrf G C). The infimum (or greatest lower bound) is 
defined in a similar manner. 

Since every element of C 2 has a corresponding supremum and infimum, the 
>Z relation defines a lattice structure on C (and on V). 




Fig. 4. The ^ relation: example. 



Figure @] gives an example of incomparable first contact 
distributions. 

It appears difficult to define a total order on C: comparing 
the distributions d\ and d 2 in figure [4] is a matter of choice 
and depends on the bundle delivery guarantees one wants to 
enforce. The > relation is thus a least common denominator, 
and could be replaced in what follows with a more restrictive 
order definition. 

The worst (smallest) element of C is the _!_ (bottom) distri- 
bution: L(t) = (\/t > 0). The best (greatest) first contact 
distribution is denoted T (top): T(i) = S(t) (Vi > 0); the S 
symbol denotes the Dirac distribution. 

The y relation can be extended to T>. For all Di, D 2 E V: 

DihD 2 <S=> VT > : D X (T, •) >r D 2 (T, •) 

The D± delivery distribution is such that VT > : 
D±(T, ■) = J_. The definition of £>t follows immediately. 

IV. Delivery distribution operators 

A. The forwarding operator 

Let D s bd be the delivery distribution associated with the 
delivery of a bundle from a source node s to a destination d 
via node b. More precisely, if s decides to send a bundle at 
time T, it will reach d after a delay described by the D s bd(T, •) 
distribution. D s m can be computed thanks to D s b and Dm'- 

D sbd = D sb <g> D bd (4) 

The <S> (or forwarding) operator is a function defined for all 
distribution pair. We have ® : V 2 — * T>: 

{D 1 ®D 2 ){T,t)^ f D 1 (T,x)D 2 (T + x 1 t-x)dx (5) 
Jo 

It is easy to see that this operator is associative but not 
commutative. 



Equation <(3j simply states that since the total delivery delay 
is equal to t, if the delay to reach b is equal to x, then the 
delay from b to d is t — x. 

Equation (O can be generalized: a bundle could be for- 
warded through several intermediate hops before reaching its 
destination. We denote D s -d (notice the dash) the delivery 
delay distribution for a bundle sent from a source s to a 
destination d at time T; from now on, <g> will thus be applied 
to any kind of delivery distributions. 

For example, the graph below depicts a simple delivery 
path, i.e. a sequence of forwarding nodes; the corresponding 
delivery distribution is also given. 

s — >- a — b — ^ d : D s - d = D sa ® D ab (g) D bd 

We say that two delivery paths with a common source s 
and destination d are disjoint if the intersection of the set of 
nodes they involve is {s,d}. 

B. The duplication operator 

Let D s ^d be the delivery distribution associated with the 

d 

delivery of a bundle from s to d if it is duplicated so as to 
follow the disjoint delivery paths described by the distributions 
D s -d and D' s d . We have: 

D*-d =D s _ d ®D' s _ d (6) 

a 

The © (or duplication) operator is a function © : T> 2 — > T>, 
defined as follows: 

(Dx®D 2 )(T,t) = (l- J D 1 (T,x)dx^J D 2 (T,t) + 

(l -J D 2 (T,x)dx\ D x {T,t) (7) 

The expected delay computed is that of the first bundle to 
reach the destination d. It is easy to see that © is associative 



and commutative. We decide that has a higher precedence 
than ©. 

Equation (0 is the sum of two terms. Each term is the 
probability that the bundle reaches the destination after a delay 
t using one path and that the bundle following the other path 
is not arrived yet. 

Notice that we have both Di®D 2 h Di and Di®D 2 h D 2 
(appendix, corollary [TJ. This means that, contrary to what 
happens in deterministic networks, duplicating a bundle to 
send it along two paths can improve performance: it is not 
the case that the best path always delivers the bundle first. 

The definition of this operator allows us to apply it to 
arbitrary independent distributions (for example, involving 
duplication and forwarding). This allows the computation of 
the distribution associated with a non trivial way to deliver 
a bundle, such as the one depicted below; the corresponding 
distribution formula is given on the right. Two arrows leaving 
a node depict a duplication. 



c ■ 



e ■ 



D s - d = (D sc D cd ) © (D sb 
(D be ®D ed ®D bf ®D }d )) 



Figure [6] shows examples of the distributions obtained 
using those operators. As expected, the "duplication" operator 
shortens the delays and increases the delivery probability. 

C. The scheduling operator 

Let D s ~d be the delivery distribution that, every time a 

'd 

bundle has to be sent, chooses the best delivery strategy out 
of D s -d and D' s d . We have: 



Ds-d =D s _ d 0D' 



S — d 



(8) 



The definition of is straightforward. It is a function 
: V 2 -> V such that: 



(Di D 2 ) (T, t) 



Di(T,t) if D 2 (T, •) ^ Di(T, •) 
D 2 (T,t) otherwise 



(9) 



If s sends a bundle at time T, it is delivered using D 2 (T, •) 
if and only if D 2 (T,-) h L>i(T, •). This operator is not 
commutative since >: is not a total order: when Z?i(T, •) and 
D 2 (T, •) cannot be compared, Di(T, •) is chosen. We decide 
that has a lower precedence than both and ©. 

The following example involves all the operators defined 
above. Two arrows leaving a node, one of them dotted, depict 
a scheduling operation. The dotted arrow leads to the second 
argument of 0, emphasizing the operator's non-commutativity. 



•/. — ./• 



D s -d = (D sc D cd ) © (D sb 

{D bf ®D }d <Z)D be ®D ed )) 



D. Delivery schemes 

We have defined a delivery path as a delivery strategy that 
only involves forwarding. 



A delivery scheme with source s and destination d is a 
general delivery strategy that allows a bundle to be delivered 
from s to d. It can use an arbitrary number of forwarding, 
duplication and scheduling operations. A delivery path is thus 
a particular delivery scheme. 

Two delivery schemes from s to d are disjoint if the 
intersection of the set of nodes they involve is {s,d}. 

V. Delivery guarantees 

Knowing the delay distribution d s - d e C associated with 
the delivery of a bundle allows us to verify a large range of 
conditions on permissible delays or on delivery probabilities. 

For example, the condition 



d s - d (t)tdt < d n 



imposes a maximum expected delay d max , while 

/ d s - d (t)dt>.9 and / d s ^ d (t) dt > .99 
Jo Jo 

matches distributions delivering a bundle in less than one hour 
nine times out of ten, and in less than a day with a probability 
of 99%. 

We naturally impose that a condition fulfilled for a certain 
delivery scheme must be fulfilled for better schemes. 

Definition 4: A delivery condition C is a predicate: 
C :C — > {true, false} with \/di,d 2 € C such that d% >z d 2 : 
C(d 2 ) C{di). 

A condition C can be extended to a delivery distribution 
D eV: C{D) VT > : C(D(T, ■)). 

VI. Delivering bundles with guarantees 

A. Probabilistic Bellman-Ford 

Algorithm Q] adapts the Bellman-Ford algorithm to pre- 
dictable disruption tolerant networks. In this section, we do not 
allow bundle duplication. Notice that, in general, the concept 
of "shortest path" is meaningless since the ^ relation is a 
partial order. 

Similarly to the Bellman-Ford algorithm, algorithm Q] com- 
putes, for every node n 6 J\f, the best distribution leading 
to the destinations found so far (B n ). This distribution is 
propagated to its neighbours (i.e. all the other nodes since 
the network is infrastructureless). 

Once node x receives the best delivery distribution B y 
found by y, it computes the delivery distribution obtained if it 
would send the bundle directly to y, and if y would forward it 
according to B y . The resulting distribution is denoted D xy - d 
(line©. 

D xy -d is compared to the best known distribution to the 
destination (B x ) by means of the operator. If D xy ^ d is 
better than B x on some time intervals, B x is updated (line [9]). 

The algorithm terminates once no more B x distribution is 
updated. 

Figures [7] and [8] demonstrate how the algorithm works by 
means of a small example. 




D (%/day) 





Fig. 6. Forwarding ((g)) and duplication (©) operators: example. We denote Di the delivery distribution depicted in figure [5] The top part of this figure 
depicts a contact profile (left) and the associated delivery distribution D2 (dark squares represent a probability equal to 1). The left 3D plot depicts D\ ® D2, 
the right one Di © D2. 



Algorithm 1: Probabilistic Bellman-Ford 



Data: d is the destination node 



1 V x G M \ {d} : B x 
i B d «- D T ; 



repeat 

stabilized <— true; 
forall x G Af do 
forall y G N do 



D 



3 
4 
5 
6 
7 
8 
9 

10 
11 

12 
13 

14 until stabilized 



D xy <g> B v ; 



xy — d 

B x ^ B x D xy - d then 
stabilized <— false; 



B x ^ B x D xy -d, 



end 



end 



end 



As mentioned before, this algorithm generalizes both [12] 
(i.e. converges to the "shortest expected path") and [170 {i.e. 
finds the exact shortest path in the case of perfectly predictable 
networks). 

The delivery computed by this algorithm depends on the 
order at which the elements of Af are picked up (lines [5] and 

4 To be fair, this work also deals with message transmission delays, which 
are not considered here. 



|6j. In practice, it might be preferable to rely on a heuristic to 
choose the preferred elements first. 

B. Guarantees 

Our aim is now to find a way to deliver bundles that 
fulfills a given condition C as specified in definition |H while 
trying to minimize the network's bandwidth/energy/memory 
consumption. 

Ideally, the DTN is predictable enough to enforce condi- 
tion C without duplicating any bundle. We thus propose to 
rely on algorithm Q] to find a first delivery scheme (and, thus, 
a first delivery distribution D\). 

If C is not fulfilled by D\, we search for another fast bundle 
forwarding scheme using algorithm [T] let D2 be its delivery 
distribution. We then duplicate the bundle on both delivery 
schemes, yielding a distribution D\ © £>2- We have already 
pointed out that D\ © D 2 h D x , thus C{D X © D 2 ) is more 
likely to be true then C{Di). 

This process is iterated until C is finally fulfilled. 

As mentioned in section IIV-BI the distribution computed 
by the "duplication" (ffi) operator is biased if its operands are 
not independent distributions. The simple distribution formula 
(D s b <g) Dm) © (D s b (g) Dbd) brings to light the problem caused 
by dependent distributions. 

To avoid this bias, we ensure that D\ and D 2 are indepen- 
dent by forbidding D 2 to rely on the nodes involved in D\ 
(source and destination nodes excluded, line |5}. 

The resulting algorithm is given below. 




This graph gives the contact profiles of the nodes 

a, b, c,d,e G TV. 

Unconnected nodes never meet each other: they have a null 
contact profile (and a corresponding delivery distribution D±). 
The label connecting the other nodes describes which days they 
might have a contact. For example, there is one chance out of 
four that b and c meet on Thursday, and one out of two on 
Friday. 



Fig. 7. A predictable network. 



B„ = Bh = B r = B,l 



D i Bp 



D 1 



The opposite table shows how our probabilistic Bellman- 
Ford algorithm behaves. This example is based on a 
simple network made of 5 nodes. The nodes contact 
profiles are given in figure [7] In this example, a is the 
source node and e is the destination. 
At first, all the nodes (but the destination) have no 
knowledge of any path to the destination; their best 
distribution is thus set to D±. The destination's delivery 
distribution to itself is of course Dt- 
Line 2 depicts the results obtained after the first iteration. 
Since only c and d have contacts with the destination, 
only B c and Bd are modified. They are set to the direct 
contact with the destination distribution since, for exam- 
ple, -D c _ e (g) D-f =D c - e . The delivery distributions are 
depicted as a square plot; the a>axis is the bundle sending 
time, the j/-axis is the delay to reach the destination. 
Each square represents a 24 hour period, the first column 
matches bundles sent on Monday. 

During the next iteration (line 3), a discovers it might 
meet with d before d meets e. B a is thus changed to 
D a _d Dd-e- The bundles received by b can be for- 
warded to c or d. The distributions D;,_ c _ e and Db-d-e 
are thus compared; bundles sent Tuesday or before are 
sent via d, those sent after Tuesday are sent via c. 
The last iteration allows a to decide when bundles should 
be sent to b or d. The distributions B a and _D a _b®_Bf, are 
thus compared; the latter is given between parentheses. 
Neither c nor d should forward bundles to b, thus B c and 
Bd are left untouched. 

The algorithm is stabilized since neither b, c, or d should 
forward bundles via a. 




B d = D d - 




B a = D a . 
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Fig. 8. Probabilistic Bellman-Ford: example. 



Algorithm 2: Constrained probabilistic delivery 
Data: Network nodes A/"; delivery condition C 
Data: Bundle source s and destination d 

1 B <- L> ± 

2 repeat 

3 
4 
5 



Using nodes in A/", compute D E V via algorithm Q] 

B ^B®D 

J\f <— N \ {nodes involved in D} U {s, d} 



6 until (7(B) or A" = {s,<i} 



Nothing guarantees of course that there exists a way to 
deliver bundles that satisfies C: even an epidemic broadcasting 
might not suffice. 



C. More on disjoint delivery schemes 

The constrained probabilistic delivery algorithm above 
computes a delivery scheme that consists of duplicating 
the bundle to multiple, independent, non-duplicating delivery 
schemes. 

To ensure independence, algorithm [2] enforces those non- 
duplicating delivery schemes to operate on completely distinct 
node sets. This might be too stringent if the network is small 
or sparse. We thus propose to allow such a delivery scheme 
to use nodes that are unlikely to receive a bundle according 
to the other schemes. The resulting delivery distributions will 
thus be almost independent. 

Line [5] of algorithm [2] is thus changed: only the nodes 
involved in D with a probability higher than a given threshold 
are removed. The specific value of this threshold is a function 
of the network considered. 

The rest of this section explains how to compute the 
probability that a given node receives a bundle, given a (non- 
duplicating) delivery scheme computed by algorithm Q] 

We have seen that the proposed modified Bellman-Ford 
algorithm does not lead to a simple routing table: if a bundle 
reaches a given node at time T, its next hop depends on its 
destination and on T. Each node n divides time in intervals 
J™, I<2, ■■ ■ (by means of the operator, algorithm Q] line |7|, 
and each interval / matches a given next hop H n (I). In the 
example figure [8] node b has defined two intervals: l\ = 
[Monday, Tuesday] and l\ — [Wednesday, Sunday]; Hb(l\) = 
d and H b {l\) = c. 

A bundle crosses a number of nodes on its way to its 
destination. We compute the probability P n that a given node 
n is one of them. 

Let s be the bundle source node and d the destination. If 
the bundle is ready to be sent at time T, it should reach n = 
H S (I), where / is the time interval of s such that T E I. 
The bundle arrival time at n follows the contact distribution 
N(x) =D sn (T,x), thus P n = J °° N{x)dx. 

Once the bundle has been received by n, each time interval 
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Fig. 9. Finding a fast delivery scheme that fulfills condition C using 
algorithms [T] and [2] This figure represents the delivery distribution lattice 
(introduced in figure [4); it is depicted the usual way (basically, an element is 
greater than another one if it is placed above). The greyed area corresponds to 
elements that satisfy condition C. 1. The adapted Bellman-Ford algorithm is 
used to find a distribution (D\ ) that characterizes a fast way to deliver bundles; 
B l s denotes the source node's best distribution found after i iterations. We have 
_!_ < Bl < Bl < ■■■ < D\. 2. Since -iC(Di), another disjoint delivery 
distribution, D2, is computed using algorithm [T] Combined with D\, it leads 
to Di©D2 which still does not satisfy C. D3 is thus computed, and combined 
with D\ and D2, gives a satisfactory delivery scheme D\ © D2 © -D3. We 
have Di H D 1 © D 2 d £>l © D 2 © D 3 . 



If matches a potential next hop ni = H n (If). We have: 
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N?(x) dx 



N(t)D nni (t,x-t)dt (10) 



(11) 



Equation ( fTOb gives the bundle time arrival distribution at 
the next hop n,. This process can be continued recursively 
until the probability of receiving the bundle is known for all 
nodes. 

The bundle forwarding process can be represented by a 
graph. The children of a node are the potential next hops. 
The graph obtained for the example depicted in figure [8] is 
given below. 
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The numbers labelling the nodes are the probabilities of 
receiving the bundle, as given by ( fTTT i. The destination e thus 
receives the bundle with a probability of ^ + jg. 

VII. Conclusion and future works 

We propose to model contacts between a disruption tolerant 
network's mobile nodes as a random process, characterized by 
contact distributions. Such a description is more general than 



those generally encountered in the literature, and allows, for 
example, to model a perfectly deterministic network. 

We show how such contact distributions can be combined to 
compute the bundle delivery delay distribution corresponding 
to a given delivery strategy (i.e. a description of the nodes for- 
warding decisions). We show how the Bellman-Ford algorithm 
can be adapted to cope with such stochastic networks. 

There is a tradeoff between a bundle's delivery proba- 
bility/delay and the consumption of network resources. We 
propose to duplicate bundles along disjoint "shortest" path so 
as to meet a given delivery guarantee without consuming too 
many resources. The corresponding algorithms are given. 

This work can be continued along several lines. 

We have proposed a way to route bundles through the 
network; other routing strategies should be explored and 
compared. 

Three operators on delivery distributions have been defined. 
Others could be added so as to describe more subtle routing 
decisions, or to deal with bundles' transmission delays. 

Real network traces should be analysed so as to quantify 
their predictability, to compare delivery strategies, and to 
measure how predictability impacts performance. 

Appendix 

Lemma 1: VDi,D 2 ,D 3 E V, we have D 2 h D 3 => D\ ffi 
D 2 h £>i ffi D 3 . 

Proof: Given the definition of ffi and y, one must prove 
that, VDi, D 2 e V, VT > 0, * > 0, given D 2 t D 3 : 
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D 2 (T,y)dy 



D\ (T, x) dx 



Di(T,y)dy 



D 3 (T,x)- 



D 3 (T,y)dy 



D 1 (T,x)dx (12) 



The left-hand part can be written as: 

/ D 1 {T,x)dx+ [ D 2 {T,x)dx- 
Jo Jo 

( ( [D 1 (T,x)D 2 (T,y) + D 1 (T,y)D 2 (T,x)]dydx 
Jo Jo 

(13) 

Changing the double integral's integration order, the last 
term of ( fT3l ) is equal to: 

/ Di(T,x) / D 2 (T,y)dydx + 
Jo Jo 

I Dx{T,x) I D 2 (T,y)dydx 

Jo Jx 

= f D 1 (T,x)dx f D 2 (T,y)dy (14) 
Jo Jo 



The same procedure can be applied to the right-hand part 
of ( fT~2l >. ( fT2] i is thus equivalent to: 

/ D 2 {T,x)dx> [ D 3 (T,x)dx (15) 
Jo Jo 

Which holds by hypothesis. ■ 

Corollary 1: VDi, -D 2 e T>, we have D x © D 2 h D x . 
Proof: From Lemma [Q D x ffi D 2 h D x © J_ = D x . ■ 

References 

[1] Z. Zhang, "Routing in intermittently connected mobile ad hoc networks 
and delay toerant networks: overview and challenges," in IEEE Com- 
munications Surveys and Tutorials, vol. 8-1, 2006, pp. 24-37. 

[2] S. Burleigh, A. Hooke, L. Torgerson, K. Fall, V. Cerf, B. Durst, K. Scott, 
and H. Weiss, "Delay-tolerant networking - an approach to interplanetary 
internet," IEEE Communications Magazine, vol. 41, no. 6, pp. 128-136, 
2003. 

[3] A. Vahdat and D. Becker, "Epidemic routing for partially connected ad 
hoc networks," Duke University, Tech. Rep. TR CS-200006, April 2000. 

[4] T. Small and Z. J. Haas, "The shared wireless infostation model: a new 
ad hoc networking paradigm (or where there is a whale, there is a way)," 
in MobiHoc '03: Proceedings of the 4th ACM international symposium 
on Mobile ad hoc networking & computing. New York, NY, USA: 
ACM Press, 2003, pp. 233-244. 

[5] P. Juang, H. Oki, Y. Wang, M. Martonosi, L. Peh, and D. Rubenstein, 
"Energy-efficient computing for wildlife tracking: Design tradeoffs and 
early experiences with zebranet," in ASPLOS, San Jose, CA, Oct. 2002. 

[6] F. Tchakountio and R. Ramanathan, "Tracking highly mobile endpoints," 
in WOWMOM '01: Proceedings of the 4th ACM international workshop 
on Wireless mobile multimedia. New York, NY, USA: ACM Press, 
2001, pp. 83-94. 

[7] A. Spuropoulos, K. Psounis, and C. Raghavendra, "Single-copy routing 
in intermittently connected mobile networks," in Proceedings of IEEE 
SECON, October 2004. 
[8] T. Spyropoulos, K. Psounis, and C. Raghavendra, "Spray and wait: An 
efficient routing scheme for intermittently connected mobile networks," 
in Proceedings of SIGCOMM 2005, 2005. 
[9] D. Nain, N. Petigara, and H. Balakrishnan, "Integrated Routing and 
Storage for Messaging Applications in Mobile Ad Hoc Networks," in 
WiOpt '03: Modeling and Optimization in Mobile, Ad Hoc and Wireless 
Networks, Sophia-Antipolis, France, March 2003. 

[10] Y. Wang, S. Jain, M. Martonosi, and K. Fall, "Erasure-coding based 
routing for opportunistic networks," in WDTN '05: Proceeding of the 
2005 ACM SIGCOMM workshop on Delay-tolerant networking. New 
York, NY, USA: ACM Press, 2005, pp. 229-236. 

[11] A. Lindgren, A. Doria, and O. Schelen, "Probabilistic routing in inter- 
mittently connected networks," SIGMOBILE Mob. Comput. Commun. 
Rev., vol. 7, no. 3, pp. 19-20, 2003. 

[12] K. Tan, Q. Zhang, and W. Zhu, "Shortest path routing in partially 
connected Ad Hoc networks," in Proc. of IEEE GLOBECOM'03, vol. 2, 
December 2003, pp. 1038-1042. 

[13] E. P. C. Jones, L. Li, and P. A. S. Ward, "Practical routing in delay- 
tolerant networks," in WDTN '05: Proceeding of the 2005 ACM SIG- 
COMM workshop on Delay-tolerant networking. New York, NY, USA: 
ACM Press, 2005, pp. 237-243. 

[14] J. Leguay, T. Friedman, and V. Conan, "DTN routing in a mobility 
pattern space," in WDTN '05: Proceeding of the 2005 ACM SIGCOMM 
workshop on Delay-tolerant networking. New York, NY, USA: ACM 
Press, 2005, pp. 276-283. 

[15] C. Shen, G. Borkar, S. Rajagopalan, and C. Jaikaeo, "Interrogation-based 
relay routing for ad hoc satellite networks," in IEEE Globecom, Taipei, 
Taiwan, November 17-21 2002. 

[16] M. Musolesi, S. Hailes, and C. Mascolo, "Adaptive routing for intermit- 
tently connected mobile ad hoc networks," in WOWMOM '05: Proceed- 
ings of the Sixth IEEE International Symposium on a World of Wireless 
Mobile and Multimedia Networks (WoWMoM'05). Washington, DC, 
USA: IEEE Computer Society, 2005, pp. 183-189. 

[17] S. Meragu, M. Ammar, and E. Zegura, "Routing in space and time in 
networks with predictable mobility," Georgia Institute of Technology, 
Tech. Rep. GIT-CC-04-7, 2004. 



